THz emission from ultrafast optical orientation 
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We show both theoreticaUy and experimentally that the magnetization density accompanying 
ultrafast excitation of a semiconductor with circular polarized light varies rapidly enough to produce 
a detectable THz field. 



A host of recent fundamental studies on the behavior of 
spins in semiconductors have been motivated by the idea 
of using the spin degree of freedom in information pro- 
cessing, the goal of a field that has come to be known as 
spintronics. Central to many of these investigations has 
been the optical injection of spin polarized electrons and 
holes into the band structure of a semiconductor. These 
spin polarized carriers are then available to be dragged 
by a bias voltage, for example, to regions of interest in 
any proposed device structure. Such optical injection 
of spin polarized carriers, also called optical orientation, 
can be effected by the absorption of circularly polarized 
light P, 0] , which takes advantage of the selection rules 
to populate particular angular momentum states of the 
band structure. 

In the analysis of most of these experiments it is as- 
sumed that the holes can be neglected, because in bulk 
semiconductors the hole spin lifetimes are typically much 
shorter than the electron spin lifetimes. This makes 
it difficult, if not impossible, to measure the optical 
orientation of hole spins in bulk semiconductors using 
usual techniques employed to study electron spin orien- 
tation. However, Hilton and Tang recently showed 
that pump-probe techniques could be used to probe the 
hole spin dynamics in bulk GaAs. They found a time of 
T/i = 110 fs for the hole spin relaxation time, which has 
spurred theoretical investigations [H . All such techniques 
measure optical orientation indirectly by monitoring, af- 
ter excitation, the population in various bands, or its 
effect on luminescence or Faraday rotation [H, [gl . 

In this manuscript we present an approach for the di- 
rect measurement and study of ultrafast optical orienta- 
tion. While we focus on GaAs, the technique is applicable 
to a wide range of semiconductors. It relies on the mag- 
netization that accompanies the injected spin density in 
optical orientation: Under ultrafast excitation this mag- 
netization varies rapidly over a sub-picosecond timescale, 
and radiates in the terahertz (THz) regime. Such ra- 
diation has already been used to study a host of other 
ultrafast processes 0, H, @, [13] j and we show here that 
detailed measurement and analysis of the emitted THz 
electric field trace can be used to study the electron and 
hole spin dynamics. 

To estimate the THz field strength radiated from ultra- 



fast optical orientation, we compare the injected magne- 
tization source M{t) in optical orientation to the injected 
polarization source P{t) in the shift current (also called 
the intrinsic photovoltaic effect [IH). When a semicon- 
ductor is excited by photons with energies above the band 
gap, electrons and holes are injected at a rate that closely 
follows the temporal pulse profile X(r, i), which is the in- 
tensity profile of the laser pulse in the sample. The car- 
rier density injection rate n is given by n(r, t) = £^i^X(t, t), 
where to refers to the carrier frequency of the laser, and 
is proportional to the absorption coefficient. The spa- 
tial coordinate r accounts for the variation from the ex- 
ponential decay into the surface, and the variation in 
intensity across the spot-size. The polarization and mag- 
netization sources we consider both arise from the carrier 
injection, although they have different symmetry charac- 
teristics and underlying microscopic origins. 

In III-V semiconductors, the shift current appears for 
excitation by certain linear polarizations of the field [l3 | . 
In this effect the center of charge of the carriers moves as 
the carriers are promoted from the valence band to the 
conduction band, and a net current directly results. 



Pinj(<) =pn{r,t), 



(1) 



where p = ed is the average injected dipole moment per 
carrier. Here, e is the electron charge and the displace- 
ment d is on the order of a Bohr radius ae. Detailed 
calculations show that for GaAs under excitation with 
linear polarized light along [110] this distance d is very 
close to the GaAs bond length, d — 2.54 ob [i3]- Since 
the polarization Pinjit) rises on a subpicosecond time- 



scale, it radiates in the THz regime [ij]. However, the 
polarization survives for much longer than a picosecond, 
and so its decay is not expected to radiate in the THz 
regime. Thus for the polarization that contributes to the 
far field radiation we have P{t) = Pinj{t). 

Optical orientation arises from circularly polarized ex- 
citation of carriers near the band edge. Because of spin- 
orbit splitting these carriers can have a net spin polar- 
ization, and accompanying this spin density is a magne- 
tization density M{r,t). For optically injected electrons, 
close to the conduction band edge, the injected magneti- 
zation density is given by 



Me,inj(r,t) = ^en(r, i). 



(2) 
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where p,e is the mean magnetic moment of an injected 
electron. Unhke for free electrons - where the magnetic 
moment of an electron is simply /i — g/j^s/h, where /^b 
is the Bohr magneton, s — h/2, and g is the free electron 
(7- factor, g — 2.0023 - the average magnetic moment of 
an injected electrons can be written as fie — glfJ-BS/h, 
where g* is an effective gr-factor of an electron, and S is 
the injected degree of spin polarization in units of h/2. 
The effective g-factor can be calculated from Roth's ex- 
pression for the magnetic moment of a Bloch state in the 
n-th band 151, 
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and defining the effective g-factor via g* = — 2/i„//XB. 
In Eq. Q , cr^ is the expectation value of the Pauli spin 
matrix in the n-th band, p^^, is the x component of the 
momentum matrix element between the bands n and s, 
and e„ is the energy of the n-th band. The z direction is 
the direction along which the magnetization is aligned. 
The appearance of cr^ in the first term is required to gen- 
eralize the more common g-factor expression (16| beyond 
j = 1/2 states. For optically injected electrons in GaAs, 
close to the conduction band edge, g^ is —0.44 at low 
temperatures fl7| . 

This magnetization density Me(t) will vary on an ul- 
trafast time-scale with the laser pulse, and it will radiate 
in the THz regime. The radiation will be magnetic-dipole 
like, in contrast to the electric-dipole like radiation from 
the shift current. Since the lifetime of the injected spins 
is known to be on the order of many picoseconds 
their decay will not contribute to the THz and we have 

Me(r,t) =Me,inj(r,t). 

The far-field radiation from both magnetic and polar- 
ization sources is dependent on the second derivatives of 
the sources, M{t) and P{t). A simple estimate of the ra- 
diation from just the electron magnetization Me{t) can 
be made by comparing this source to the polarization in- 
duced by the shift current. Since both of these sources 
survive long enough that their decay profiles are not ex- 
pected to radiate in the THz regime, it is mainly the 
generation (or turn on) of these sources that provides the 
THz radiation. The ratio of the peak strengths for the 
magnetization and polarization sources are directly com- 
parable, since the different Fresnel coefficients governing 
the transmission of radiation for each type of source only 
determine the radiation patterns. So to first order we 
can estimate the strength of the magnetic radiation from 
the ratio Me{t)/P{t). Using Eq. dH) and ^ gives 
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where a is the fine structure constant. In GaAs optical 
orientation gives 50% spin polarization [S = 0.5?i/2), so 
the ratio Me/P is roughly 1.6 x lO"*. The THz fields 



from the shift current have been measured to be about 
1500 V/m, so an absolute sensitivity of about 0.2 V/m 
is needed to measure the THz radiation from optical ori- 
entation of bulk GaAs electrons alone. This requirement 
is experimentally challenging |il8i] but would not be im- 
possible. Note that this is because the THz field ampli- 
tude, rather than the intensity, is measured experimen- 
tally. Were the intensity measured, the signal from the 
magnetization effect would be down from that of the shift 
current by a factor of {Mf,/PY ■ 

The relatively short lifetime of the holes, however, sug- 
gests that the hole spins may radiate in the THz regime. 
Including the decay of the hole spins, the magnetization 
of the holes follows 



Mft(r,t) = Af,,,inj(r,t) 



where M^^inj is the hole injection rate 



(5) 
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analogous to the electron injection rate, but here flh is 
the mean magnetic moment of the injected holes. 

On the one hand, the decay of the holes adds more 
temporal variation to the magnetization, and enhances 
the THz signal; on the other hand, since the decay lowers 
the overall magnitude of the magnetization, it is expected 
to decrease the THz signal. To include both effects, we 
have solved for the far-field radiation from an injected 
magnetization using a Green function formalism tailored 
to planar interfaces [l^. We find that for pulse widths 
of 100 fs the effect of the hole spin decay {rh = 110 fs) 
is to decrease the peak THz radiation from this source 
by roughly a factor of two. That is, the hole spins decay 
fast enough so that the decay process will contribute the 
THz radiation, but slow enough that the overall signal 
from it is not washed out. 

Yet the relatively short lifetimes of holes, in bulk mate- 
rials such as GaAs, prevents the study of their dynamics 
by many other techniques. Thus the hole g-f actors, and 
their magnetic moments in general, have been largely ne- 
glected in bulk materials. In GaAs based nanostructures 
however, hole g-f actors have been the subject of much 
recent attention [2^. It is typically found that in those 
systems that the holes have a g-factor many times larger 
than the electron (7-f actor. In this work however, we are 
interested in the optically injected magnetic moments 
/i/i, rather than the g-factors themselves. The average 
magnetic moment flh will depend on the details of the 
populations injected into the heavy-hole and light-hole 
bands. We calculate directly the injected hole magneti- 
zation into these states, but for easy comparison to the 
electron magnetization, we quote our result as a hole g- 
factor 5^ relative to the injected electron spin density, so 
that 
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FIG. 1: (color online) (a) Experimental setup, (b) Top view 
of sample with epoxy hemisphere. PEM: photoelastic mod- 
ulator, HWP: half waveplate, QWP: quarter waveplate, PD; 
photodiode 



The tilde is to emphasize that we are not defining the 
g-factor relative to the hole spins. Following a simi- 
lar approach as in [22] we use Fermi's golden rule to- 
gether with Roth's formula Eq. ^ for the magnetization, 
and use a 30-band k • p band structure [2l| to evaluate 
Mh^inj directly. We find that 5^ — 14.1, which implies 
that the hole magnetization injection is approximately 
32 times larger than the electron magnetization injection 
and clearly dominates it [23]. 

Combining this with our calculations for the radiation 
emission (following [l^) we find that for GaAs the hole 
radiation should be of the order 3 x 10^"^, or about 1/300 
times the shift current radiation. 

Despite the relatively large expected signal from the 
the holes, the ratio of THz emission from magnetization 
to shift current is still small, and one must choose an 
experimental geometry which minimizes the shift cur- 
rent. Differences in the spatial symmetry characteristics 
of shift current and magnetization sources are exploited 
to distinguish the weak magnetization source from the 
much stronger shift current source. 

Fig. [T] illustrates the experimental setup for observ- 
ing the transient THz radiation. An 80 MHz Ti:sapphire 
oscillator delivers InJ, 100 fs pulses at 800 nm. For this 
wavelength, the injected kinetic energy of heavy (light) 
holes is 3.6 meV (9.5 meV), much less than the 300 meV 
needed to excite holes in the spin-orbit split-off band. 
The pulses are split into a pump beam for sample exci- 
tation and a probe beam to sample THz radiation via 
electro-optic sampling. The pump beam is incident onto 
a 100 /im thick (110) GaAs sample with a spot diam- 



eter of 25 /Ltm, producing a peak incident intensity of 
500MW/cm^ and a carrier density of 2 x 10^^ cm~'^ over 
an absorption depth of about 1 /im. Since the total mag- 
netization or electric polarization is determined by total 
number of carriers excited, and not their density, there is 
no premium associated with beam focusing. The emitted 
THz radiation emitted along the surface normal direction 
from the sample back side is collected with a f = 50 mm, 
90° off-axis parabolic and focused with the same type 
of mirror onto a 500 /im thick (110) ZnTe electro-optic 
sampling crystal. The 800 nm probe pulse is spatially 
overlapped with the THz field inside the ZnTe crystal 
and temporally scanned. The pump beam polarization is 
controlled by half- and quarter-wave plates (QWP) and 
a photo-elastic modulator (PEM). The half- wave plate 
(HWP) is oriented to provide light linearly polarized 45° 
from the s polarization vector (the vertical). The PEM 
operates in quarter-wave mode and modulates the polar- 
ization state at a frequency of 42 kHz, enabling lock-in 
detection techniques. 

A major difficulty in measuring the emitted THz is 
that the injected magnetic dipoles are parallel to the laser 
beam propagation direction. For normal incidence exci- 
tation, this injects magnetic dipoles that are perpendic- 
ular to the surface. The component of the emitted radi- 
ation nearly parallel to the surface is small, and so there 
is little transmission. In the shift current effect however, 
the polarization is parallel to the surface, and this allows 
the radiated THz field to easily transmit through the sur- 
face. To enhance the emission of THz from the ultrafast 
optical orientation it is thus necessary to irradiate the 
sample at a high angle of incidence; we use an angle of in- 
cidence of 55°, exciting through an epoxy hemisphere (see 
Fig. [TJd), to give an angle of refraction inside the GaAs, 
6r, of about 20°. Even for this small angle it follows 
from the Maxwell saltus (boundary) conditions that the 
emitted THz radiation mainly propagates in the direction 
of the surface normal, because the induced electron-hole 
plasma is < 1 /im thick. 

We define /3 as the angle between the QWP's fast axis 
and s. The shift current and magnetization signals both 
depend on /? in a characteristic and unique way, so we 
focus on that variable to distinguish between the two 
THz sources. If we define 7 as the angle between the 
[001] crystal axis (lying in the plane of the sample) and 
s, the internal THz field emitted by the two types of 
sources are given by 

AEf^^ = Cedcos(6'H)(2sin7-3sin3 7)cos(2/3),(8) 



C/ieffsin {Or) sin^ (2/3), 



(9) 



where the proportionality constant C includes universal 
constants and Fresnel reflection and transmission coef- 
ficients. The effective magnetic moment fieff describes 
the sum of electron and hole magnetic moments, al- 
though we expect it to be dominated by hole term. It 
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FIG. 2: (color online) (a) Simulation of THz signal vs QWP 
angle /3. red dashed: shift current signal, blue thin: mag- 
netization signal, black thick: total signal (b) Experimental 
results (points) with Fourier series fit (curve) 



is evident that the shift current (magnetization) sig- 
nal is maximized (minimized) when /3 = 0° and mini- 
mized (maximized) when [3 — 45°. This is because when 
cos(2/?) = (sin(2/3) = 0), the light polarization is mod- 
ulating between -|-45° and —45° linear (right and left cir- 
cular) states. We choose 7 = 0.02 rad (1.2°), for which 
2 sin 7 — 3sin'^7 ~ 27, in order to provide a small shift 
current source with which to compare the magnetization 
signal. 

Fig. [2^ illustrates the expected variation of the THz 
field strength from the shift and magnetization sources, 
taking the ratio of the source strengths to be 0.015. Other 
possible sources of THz radiation were neglected. Opti- 
cal rectification is small relative to the shift current effect 
above the band gap, and shares the same symmetry de- 
pendence. Sources of THz radiation related to Dember 
fields have no dependence on crystal orientation or pump 
polarization. Fig. [2}d shows results from experimental 
measurements of the THz field strength as a function of 
/3. The THz field strength is defined as the peak-to-peak 
value of the observed THz trace, although several other 
definitions of signal strength were used and each resulted 
in a similar pattern as a function of /3. 

The deviation of the expected sinusoidal pattern cre- 
ated by the shift current alone is attributed to the signal 
from the magnetization source. The zero crossings of the 
Fourier series fit to the data of Fig. [2b are separated by 
more than 90°. This behavior indicates the presence of an 
additional THz source with a /3-dependence as shown in 
Fig. [2^. The difference between maximum and minimum 



signal values recorded at /3 = 90° and 0° respectively was 
0.05 ± 0.14 /xV. This difference would have had to be on 
the order of 1 /xV to attribute the data's deviation to a 
uniform vertical shift. Experimental data indicates that 
the relative strengths of the magnetization and shift cur- 
rent sources have a ratio 0.017, which is within a factor of 
six of the theoretical value of the hole contribution. More 
work is required to confirm that the offset in Fig. [2)3 is 
mainly due to the magnetization source. Given all the 
possible sources of systematic error, we consider this ex- 
tremely good agreement in the first study of this effect. 

Increasing the THz emission from the spin injection 
can be achieved in a few ways: The THz generation de- 
pends on the injection rate of spin polarized carriers, and 
so using a shorter pulse will give a larger injection rate 
and increase the THz signal. One could also use materi- 
als with a larger 5-factor. Materials with (7- factors near 
50, such as InAs, have been reported [3]. Finally, some 
materials also allow for a higher degree of spin polar- 
ization S to be achieved; in strained GaAs or the II- VI 
wurtzitc materials, for example, optical orientation can 
give an almost 100% polarized electron spin population. 

In conclusion, we have investigated the possibility of 
measuring THz radiation from the ultrafast optical ori- 
entation of electrons and holes. This THz radiation is due 
to a rapidly-varying magnetic source, and we have shown 
that it is large enough to be detected against the more 
usual THz signal from electric current sources. Measure- 
ments of the THz radiation from spin-injection can com- 
plement the wide range of experiments on spin polariza- 
tion decay already existing in the literature, providing a 
direct measurement of the injected spin dynamics of both 
electrons and holes in semiconductors. 
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